Leptonic decay-constant ratio f K + / from lattice QCD with physical light quarks 



A. Bazavov, 1 C. Bernard, 2 '* C. DeTar, 3 J. Foley, 3 W. Freeman, 4 Steven Gottlieb, 5 
U.M. Heller, 6 J.E. Hetrick, 7 J. Kim, 8 J. Laiho, 9 - 10 L. Levkova, 3 ' 8 M. Lightman, 2 
J. Osborn, 11 S. Qiu, 3 R.L. Sugar, 12 D. Toussaint, 8 ^ R.S. Van de Water, 13 ' * and R. Zhou 5 

(MILC Collaboration) 

'Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 
Department of Physics, Washington University, St. Louis, MO 63130, USA 
3 Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84-112, USA 
^Department of Physics, The George Washington University, Washington, DC 20052, USA 
5 Department of Physics, Indiana University, Bloommgton, IN 47405, USA 
6 American Physical Society, One Research Road, Ridge, NY 11961, USA 
7 Physics Department, University of the Pacific, Stockton, CA 95211, USA 
8 Department of Physics, University of Arizona, Tucson, AZ 85721, USA 
SUPA Department of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, Scotland, UK 
10 Department of Physics, Syracuse University, Syracuse, NY 13244, USA 
11 Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, IL 60439, USA 
Department of Physics, University of California, Santa Barbara, CA 93106, USA 
13 Theoretical Physics Department, Fermi National Accelerator Laboratory, Batavia 60510, USA 

(Dated: January 28, 2013) 

A calculation of the ratio of leptonic decay constants f K + / f 7T + makes possible a precise deter- 
mination of the ratio of CKM matrix elements IK^I/IKidl in the Standard Model, and places a 
stringent constraint on the scale of new physics that would lead to deviations from unitarity in the 
first row of the CKM matrix. We compute f K + / f n + numerically in unquenched lattice QCD using 
gauge-field ensembles recently generated that include four flavors of dynamical quarks: up, down, 
strange, and charm. We analyze data at four lattice spacings a « 0.06,0.09,0.12, and 0.15 fm with 
simulated pion masses down to the physical value 135 MeV. We obtain f K + / f 7T + = 1.1970(29)(57), 
where the errors are statistical and total systematic, respectively. This is the first physics result from 
our Nf = 2 + 1 + 1 ensembles, and the first calculation of f K + / f^+ from lattice-QCD simulations at 
the physical point. Our result is the most precise lattice-QCD determination of f K +/ 'f w +, with an 
error comparable to the current world average. When combined with experimental measurements of 
the leptonic branching fractions, it leads to a precise determination of |K,a|/|Viid| = 0.2304(12)(4) 
where the errors are theoretical and experimental, respectively. 

PACS numbers: 12.38. Gc, 12.15.Hh, 13.20.-v 



Motivation. - Leptonic decays of charged pseu- 
doscalar mesons are sensitive probes of quark flavor- 
changing interactions. Experimental measurements of 
the leptonic decay widths, when combined with pre- 
cise theoretical calculations of the leptonic decay con- 
stants, enable the determination of elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix. Further, when combined with independent determi- 
nations of the CKM matrix elements from other processes 
such as semileptonic meson decays, they make possible 
precise tests of the Standard-Model CKM framework. 
Here we present a lattice-QCD calculation of the decay- 
constant ratio fx+ 1 fir+ , which may be used to determine 
\V us \/\V ud \ via [1, 2] 

TjK -> in) |K, I 2 P K vfiK j 1 zMj) n | x i 
r(7r -+ i n ) |Kdl 2 / 2 ^ _ r£ y [i + ° EMl ' 

(1) 

where i = e, fj, and Sem is a known 1-loop QED cor- 
rection [3, 4]. The ratio fK+lf-K+ can be calculated to 
subpercent precision using numerical lattice-QCD simu- 



lations [5, 6] because the Monte Carlo statistical errors 
are correlated between the numerator and denominator 
and the explicit dependence on the lattice scale drops 
out. 

Lattice QCD enables ab initio calculations of non- 
perturbative hadronic matrix elements. The QCD La- 
grangian is formulated in discrete Euclidean spacetime, 
after which the now finite-dimensional path integral can 
be solved numerically with Monte Carlo methods. Once 
the Nf + 1 parameters of the QCD action (quark masses 
and gauge coupling) are fixed by matching to an equal 
number of experimental measurements, all other out- 
puts of the lattice simulation are predictions of the the- 
ory. In practice, currently available computing resources 
limit the lattice spacing in most recent simulations to 
a > 0.06 fm, the spacetime volume of the lattice to 
L 3 ~ (4 fm) 3 , and the pion mass to > 200 MeV. 
These choices introduce systematic errors that must be 
quantified, but the good agreement between lattice-QCD 
calculations and experiment for a wide range of observ- 
ables [7, 8] demonstrates that the errors are controlled. 

Many precise lattice-QCD weak-matrix element calcu- 
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lations relevant for quark-flavor physics have been ob- 
tained using our Nf = 2 + 1 flavor gauge-held ensembles 
generated using the asqtad-improved staggered quark ac- 
tion for the dynamical u, g?, and s quarks [9—11]. Our 
calculations of Jk+I ' fn+ on those ensembles [2, 11, 12] 
ultimately reached a precision of ~ 0.6%. We are now 
embarking on a new program of configuration genera- 
tion using the highly-improved staggered quark (HISQ) 
action [13], which has two primary advantages. The 
most important staggered discretization errors are ap- 
proximately 3 times smaller for the HISQ action than for 
the asqtad action at the same lattice spacing [14]; this 
makes it possible to undertake simulations at the physical 
pion mass, which require a box of approximately (5 fm) 3 . 
The improved quark dispersion relation for the HISQ ac- 
tion also allows the inclusion of dynamical charm quarks 
in the simulations. 

One of the largest sources of uncertainty in lattice- 
QCD determinations of quantities such as fx+ is 
from the extrapolation of the simulation results to the 
physical up- and down-quark masses. In this work, we 
replace this error with a negligible interpolation error by 
simulating at light-quark masses very close to or even 
below their physical values. The calculation presented 
here, using the configurations with four flavors of dy- 
namical HISQ quarks, lays the foundation for a larger 
quark flavor-physics program to compute pion, kaon, 
and charmed weak matrix elements, and eventually even 
those with bottom quarks, on these ensembles. 

Lattice-QCD Calculation. — This calculation is based 
on a subset of our (2 +1+1) -flavor ensembles described in 
Refs. [14, 15] . These ensembles use a one-loop Symanzik- 
improved gauge action for the gluons [16, 17], and the 
HISQ action for the dynamical u, d, s, and c quarks. The 
simulated strange and charm sea-quark masses (m s and 
m c ) are fixed at approximately their physical values. The 
simulated up and down sea-quark masses are taken to be 
degenerate with a common mass mi ps m s /27, such that 
the simulated pion mass is approximately the physical 
value 135 MeV [18]. Additional ensembles with slightly 
heavier pions (mi = m s /10) allow us to correct fx+l ' fn+ 
a posteriori for the slight mistuning of the simulated sea- 
quark masses. The spatial volumes are sufficiently large 
(m^L > 3.5) that finite- volume corrections are expected 
to be at the subpercent level for light pseudoscalar meson 
masses and decay constants, and we confirm this with an 
explicit finite-volume study. We use four lattice spac- 
ings a fa 0.15, 0.12, 0.09, and 0.06 fm to enable a con- 
trolled extrapolation to the continuum (a — > 0) limit. 
The specific numerical simulation parameters are given 
in Table I. 

Naive lattice discretizations of the quark action suf- 
fer from the problem of fermion doubling. The stag- 
gered quark action only partly eliminates fermion dou- 
blers, reducing the number of species from 16 to 4. The 
remaining species are referred to as "tastes" to distin- 



TABLE I: Nf = 2 + 1 + 1 gauge-field ensembles used in this 
work. The lattice spacing in fm is approximate. Ensembles 
for which the number of equilibrated lattices JVi ats has reached 
1000 are considered complete. At a pa 0.12 fm, we have three 
ensembles with identical simulation parameters except for the 
spatial volume. For these ensembles, the given pion masses 
and A r i ats are approximate, since they differ slightly for the 
three volumes. 



a(ftn) 


L(fm) 


M XjS (MeV) Mir 


,RMs(MeV) Mats 


0.15 


3.66 


214 


352 


1000 


0.15 


4.82 


133 


311 


1000 


0.12 


{3.94, 4.96, 5.83} 


215 


294 


1000 


0.12 


5.83 


133 


241 


1000 


0.09 


4.33 


215 


244 


1000 


0.09 


5.63 


128 


173 


707 


0.06 


3.79 


223 


229 


678 


0.06 


5.45 


134 


143 


192 



guish them from physical flavors. Quarks of different 
taste can interact by exchanging gluons with momen- 
tum components close to the lattice cutoff p^ ~ n/a. 
Taste-changing interactions split the mass degeneracy of 
pions composed of different quark tastes. The taste split- 
tings are of (D(a 2 s a 2 ) for the HISQ action, and decrease 
rapidly with lattice spacing. The taste-Goldstone and 
root-mean-squared (RMS) sea-pion masses are given in 
Table I; the difference between the two is only about 
10 MeV for the finest a ss 0.06 fm ensembles used in this 
analysis. More details can be found in Ref. [15]. 

The dynamical HISQ simulations use the fourth-root 
procedure to remove the unwanted taste degrees of free- 
dom. Due to taste-changing interactions, the rooted 
Dirac operator is nonlocal at nonzero lattice spacing, 
and leads to violations of unitarity [19-22]. Nevertheless, 
there are strong theoretical arguments and numerical ev- 
idence that the continuum limit of the rooted staggered 
lattice theory is indeed QCD [23-30]. 

We obtain the meson masses and decay constants from 
fits of two-point correlation functions with a pseudoscalar 
interpolating operator at both the source and sink. For 
each valence-quark mass combination, we compute the 
decay constant using a partially-conserved axial current 
relation: 

F xy = (m x + m y )(0\x l5 y\P X y)/M^ y , (2) 

where M xy is the mass of the pseudoscalar meson P xy 
with bare valence-quark masses m x and m y . We choose 
a fitting range for the correlators such that, in the case 
of degenerate valence-quark masses, a single-exponential 
fit gives a good fit as measured by the correlated x 2 /dof 
and p-value. We include both the ground state and an 
opposite-parity excited state for the nondegenerate cor- 
relators. Statistical errors are estimated by jackknifing. 
The ground-state mass and amplitudes are stable with 
respect to a reasonable reduction of the minimum time 
separation of the source and sink, and with the inclusion 
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FIG. 1: Pion mass (left) and decay constant (right) versus lattice spatial extent. The open circles show the numerical data 
with statistical errors. The staggered NLO xPT prediction is shown as a solid red line. The continuum NLO (NNLO) xPT 
predictions are shown as dashed (dot-dashed) blue lines. The close agreement of the staggered and continuum NLO xPT fits 
for f w is a numerical coincidence resulting from the cancellation of effects from pion taste splittings and quark-disconnected 
hairpin contributions, and is not seen at other lattice spacings. 



of additional excited states; therefore the error due to 
excited-state contamination can be neglected compared 
with the statistical errors. Additional details of these fits 
can be found in Ref. [31]. 

Chiral perturbation theory (xPT) predicts that the 
finite-volume corrections to pion and kaon masses and 
decay constants are at the per mill level in our simu- 
lations, and are therefore comparable to the size of the 
statistical errors in our numerical data. We check this ex- 
pectation with an explicit comparison of these quantities 
measured on three a ~ 0.12 fm ensembles with identical 
simulation parameters except for the spatial lattice vol- 
ume (see Table I). We fit the data at the three volumes 
to the functional form: 

aX(L)=A x (l+B x (m 7r ,L)) (3) 

where X = {m^, f w ,mx, fx}, A x is a free parameter, 
and the function B x is determined at a given order in 
%PT. We try both the NLO expression for B x in stag- 
gered xPT [32, 33], and the continuum NLO and re- 
summed NNLO expressions [34]. Figure 1 compares the 
numerical lattice data for m v and f n with these expres- 
sions. NLO staggered xPT describes the f v data very 
well at all three volumes, and describes the m T , rrix, 
and fx data adequately (to ^0.4% or better). We there- 
fore use NLO staggered finite- volume xPT to correct our 
simulation data in our central fit, and use the contin- 
uum NLO and NNLO finite-volume xPT corrections to 
estimate the systematic uncertainty. 

Because our lattice-QCD simulations are isospin- 
symmetric (the up and down sea-quark masses are equal) , 
we adjust the experimental inputs to what they would 
be in a world without electromagnetism or isospin viola- 
tion before matching the simulation data to experiment 
to find the strange quark mass m s and the average light 
quark mass to — (m u +m<j)/2. We follow the approach of 
Ref. [2] , but with updated values for the EM correction to 



Dashen's theorem. For the kaon, we consider the isospin- 
averaged mass Mi = (M% + + M|. )qcd/2, where the 
subscript "QCD" indicates that the leading EM effects in 
the masses were removed from the experimental masses 
[18]; we take the parameter A EM = 0.65(7) sta t(14)syst 
that characterizes violations of Dashen's theorem from 
our ongoing lattice QED+QCD simulations using asqtad 
sea quarks [35]. 

In this work, we tune the quark masses and lattice 
scale and determine //„■+ in a single, self-contained 
analysis as follows. We begin with numerical lattice data 
for meson masses and decay constants for a selection of 
valence-quark masses that allow us to adjust for slight 
mistunings and interpolate or extrapolate to the physical 
values. In order to reduce autocorrelations below mea- 
surable levels, the single-elimination jackknife results are 
blocked by 20 lattices (100 or 120 molecular dynamics 
time units) before the final averaging. The physical- 
quark-mass a « 0.06 fm ensemble has an insufficient 
number of configurations for blocks that large. Instead, 
we block by 5 trajectories on that ensemble and rescale 
the resulting errors by a factor of 1.17, which is deter- 
mined by comparing the errors at block sizes around 20 
to those at block size 5 on the to; = 0.1to s , a ~ 0.06 fm 
ensemble, as well as an additional to/ = 0.2to s , a « 0.06 
fm ensemble not otherwise discussed here. We then con- 
struct the ratio of squared pseudoscalar-meson mass to 
squared decay constant, M xx /F xx , which xPT predicts 
to be approximately linear in the valence-quark mass m x . 
On each ensemble, we interpolate or extrapolate data for 
this quantity at the three lightest valence-quark masses 
quadratically in m x to where it equals the experimen- 
tal value for M% /f%+- This fixes the physical value of 
to. We obtain the lattice spacing from requiring that 
the decay constant at this mass equal the experimental 
value for /„.+ [36]. We then take nondegenerate "kaons" 
in which the strange valence-quark mass (m y ) is 1.0 or 
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FIG. 2: Continuum extrapolation of /k+//tt+- Our central 
value is obtained from a quadratic-in-a 2 fit to the physical- 
quark-mass data points (red circles), adjusted for tuning er- 
rors to the physical sea quark masses using the 0.1m s data 
points (blue squares) . The size of the adjustments is too small 
to be visible on this plot. The black star shows fx+/fn+ 
at the physical quark masses in the continuum, where the 
inner solid error bar is statistical and the outer dotted er- 
ror bar includes the continuum-extrapolation error added in 
quadrature. An alternative, constant-in-a 2 fit to the finest 
two physical-mass data points, is shown by the dashed hor- 
izontal line. Its deviation from the central value gives our 
estimate of the continuum extrapolation error (see text). 



0.8 times the strange sea-quark mass and the lighter 
valence-quark mass (m x ) is the lightest available, and 
linearly interpolate in the valence strange-quark mass to 
where 2M 2 y - M xx = 2M| - M 2 ; this fixes the phys- 
ical strange-quark mass m s . Once we know m s , we ob- 
tain the u-d mass-splitting from the difference between 
the EM-subtracted neutral and charged kaon masses: 
m d - m u = {M\ a - M?)qcd x dM% x /dm x . Finally, we 
obtain f K + by linearly interpolating the decay constant 
F xy in the light valence-quark mass to m u and in the 
strange valence-quark mass to m s . In f K + , we neglect 
isospin violations from the sea quarks, and we neglect 
all isospin violations in /„.+ because these effects are of 
NNLO in xPT [6]. We also slightly correct fx a poste- 
riori for sea-quark mass mistuning using the slope with 
respect to the light sea-quark mass, dfx/dmi, measured 
by comparing the value of fx from the physical-mass en- 
sembles with that from ensembles with mi ss 0.1m s . The 
results for f K + //„■+ are shown versus a 2 in Fig. 2. 

We extrapolate f K + / f v + at physical quark masses to 
the continuum limit quadratically in a 2 . We estimate 
the uncertainty in the continuum extrapolation from al- 
ternative fits with different ansatze for the a 2 dependence 
and excluding the coarsest ensemble(s). Two reasonable 
choices are a linear extrapolation of the three physical- 
mass data points with a < 0.12 fm, and a constant fit to 
the two physical-mass data points with a < 0.09 fm. Of 



these, the constant fit, shown in Fig. 2, gives the larger 
difference from the central extrapolation, 0.43%. Ad- 
ditional extrapolations have also been tried, including a 
quadratic-in-a 2 , linear in m; fit to all the mi rs m s /27 and 
mi = 0.1m s data points, and a similar fit restricted to en- 
sembles with a < 0.12 fm. These fits give differences from 
the central value that are comparable to the 0.43% dif- 
ference seen above, as does a simple linear-in-a 2 extrapo- 
lation of the two finest a < 0.09 fm physical-mass points. 
We thus take 0.43% as our estimate of the discretiza- 
tion error. Note that, when the scale-setting procedure 
is changed by repeating the analysis using other quanti- 
ties to fix the lattice scale instead of (see Ref. [15] for 
further details), we find a 0.17% difference with the cen- 
tral value, well within the estimate for the continuum ex- 
trapolation error. To estimate the finite-volume error, we 
repeat the entire analysis using finite- volume corrections 
for the pion (kaon) masses and decay constants calculated 
at NNLO (NLO) in continuum xPT. We propagate the 
statistical uncertainties in the tuned quark masses and 
lattice scale throughout the analysis via the jackknife 
method, so they are already folded into the statistical 
error. The systematic uncertainties from the interpola- 
tions/extrapolations that fix m and m s to their physical 
values are negligible. We estimate the uncertainty due to 
EM effects by varying the values of the EM-subtracted 
meson masses used in the quark-mass tuning; this pri- 
marily affects m u . We vary the parameter Aem by its 
total error [35]. We also repeat the tuning including a 
previously-neglected EM correction to the neutral kaon 
mass, SM^o EM = 901(8) stat MeV 2 [37]. Note that direct 
EM effects on the weak matrix elements are by definition 
removed from f w + and fx+ [18]. 

Results and Conclusions. — We obtain the following 
determination of the ratio fx+ 1 fir+ '■ 

f K +/U+ = 1.1970(29) st at(54) a2 (17) F v(2) E M, (4) 

with an approximately 0.5% total uncertainty. Our result 
agrees with, but is more precise than, other independent 
lattice-QCD calculations [12, 38-41], and its error is com- 
petitive with that of the current lattice-QCD world aver- 
age [5, 6]. Because we use f„ to set the lattice scale, we 
also obtain a result for f K + = 156.10(38) stat (74) sys (24)^ 
where the errors are due to statistics, systematics in 
/if/Zwi and the uncertainty in f n . This agrees with the 
experimental determination of f K + assuming CKM uni- 
tarity [36]. 

By combining fx+ / f„+ from Eq. (4) with recent exper- 
imental results for the leptonic branching fractions [42], 
we obtain \V us \/\V ud \ = 0.2304(12) thco .(4) cxp .. Taking 
\V u d\ from nuclear (3 decay [43], we also obtain |V^ S | = 
0.2245(12)thoo.(4)exp.(l)y„ (J , which agrees with and has 
comparable errors to the determination from K — > ixlv 
semileptonic decay [42, 44-46]. Further, our result places 
stringent constraints on new-physics scenarios that would 
lead to deviations from unitarity in the first row of the 
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CKM matrix. Wc find 1 - \V ud \ 2 - \V US \ 2 - \ v ub \ 
—0.0004(7), which is in excellent agreement with the 
Standard-Model value of zero. 

We are currently extending the ensemble with the 
physical pion mass at a ~ 0.06 fm as well as the other en- 
sembles with fewer than 1000 equilibrated lattices. We 
will include this additional data in a longer work that 
updates this result and presents determinations of the 
charmed pseudoscalar decay constants. We will also per- 
form a more sophisticated analysis using staggered chi- 
ral perturbation theory [32, 33, 47] to obtain the low- 
energy constants of x?T. Our result for /ir+// 7r + lays 
the foundation for a new large-scale lattice-QCD physics 
program using the Nf = 2 + 1 + 1 flavor HISQ gauge-field 
ensembles with physical pion masses. These ensembles 
will enable calculations of light and heavy-light hadronic 
weak matrix elements with unprecedented precision, and 
will ultimately strengthen tests of the Standard Model 
and improve the reach of searches for new physics in the 
quark-flavor sector. 
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